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Abstract 

It is believed that most quasars and galaxies present two common features: the 
presence in their core of a supermassive object, and the experience of one or more 
encounters with other galaxies. In this scenario, it is likely that a substantial fraction 
of active galactic nuclei harbour a supermassive binary, fueled by an accretion disk. 
These binaries would certainly be among the strongest sources of sinusoidal gravitational 
waves. We investigate their evolution considering, simultaneously, the accretion of the 
black hole's masses from the disk, and the gravitational waves emitted during the orbital 
motion. We also consider other astrophysical scenarios involving a coalescing binary 
with non constant masses. 



1 Introduction 

Soon after the discovery of the astrophysical phenomena subsequently indicated with the 
name of Active Galactic Nuclei (AGN), it was argued that their power supply was ultimately 
gravitational in origin. If this is the common feature of the wide range of models included 
in the AGN category, then the natural conclusion is, as first pointed out by Zeldovich &; 
Novikov (1964) and Salpeter (1964), that an AGN 'prime mover' is a supermassive black 
hole (SBH). Among other compelling arguments in favour of the black hole hypothesis are 
their efficiency and stability, together with various observational discoveries, like rapid X- 
ray variability, small scale radio jets, broad emission lines, etc. (Blandford 1990; Osterbrock 
1993). 
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For these reasons, most of the theoretical work about the AGN phenomenon has been 
focused on the SBH hypothesis, namely that essentially all active galaxies contain ~ 10^ — 
10^ Mq black holes in their nuclei, and that these objects, together with their orbiting 
accretion disks, are the prime movers for most of the powerful activity. 

On the other hand, there are also compelling reasons to believe that a great number of 
galaxies have undergone at least one merger since the epoch of their formation (see Rees 
1990, and references therein). Indeed, many current models assume that the central object is 
activate, or simply refueled, as a result of these interactions with another galaxy (Osterbrock 
1993). This new evidence, along with the central SBH hypothesis, suggests that, due to 
the dynamical friction exerted by the surrounding environment, a certain number of active 
galaxies could harbour a binary black hole system (BBHS), with separation of the order 
of parsecs. This conclusion is supported by the observed bending and apparent precession 
of radio jets emerging from AGN (Begelman, Blandford, &; Rees 1980). In fact, the S- 
symmetry observed in many radio sources and in a considerable fraction of quasars at 
z < 1 (Hutchings, Price, & Gower 1988) might be due to the presence of such binaries. 
However, in many quasars the jets are strongly curved on the milliarcsec scale as well. If 
this curvature is also due to precession, much shorter precession periods, and thus smaller 
binary separations, are required. This indicates that mass flow into the galactic nucleus 
through the accretion disk dominates not only the activity, but also the evolution of the 
central system (Roos 1989). 

Moving from these considerations, recent works have investigated the evolution of a 
BBHS in the violently relaxed core of a merged galaxy, taking into account the flow of 
gas and stars into the newly formed nucleus (Begelman et al. 1980; Roos 1989; Ebisuzaki, 
Makino, & Okumura 1991; Fukushige, Ebisuzaki, & Makino 1992). 

These AGN models are interesting per se, and also in conjunction with the issue of 
gravitational waves (GW) detection. In fact, close binary systems of supermassivc compact 
objects are currently considered as the most certain observable sources for detectors like 
LIGO (Abramovici et al. 1992) and VIRGO (Bradaschia et al. 1990), and also for those 
experiments based on the Doppler tracking of an interplanetary spacecraft (see, e.g., Thorne 
& Braginsky 1976; Bcrtotti ct al. 1992). In the standard model for the emission of GW 
from a binary system, the energy loss is proportional to the square of the third derivative 
of the quadrupole moment. This description, however, does not take into account other 
energy loss mechanism, like friction with an accretion disk. 

In this work, we analyze the evolution of a BBHS when the effects of mass accretion and 
GW emission are considered simultaneously. In §2 we will describe in detail a single BBHS, 
and find, under some reasonable assumptions, the behavior of the separation between the 
two components. In §3 we show how the detection of these waves could provide some useful 
information on the physical characteristic of the AGN, namely its mass and accretion rate. 
In §4 we extend our analysis to a simple population of BBHS systems, and find the evolution 
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law for the distribution function. We conclude, in §5, pointing out other astrophysical 
situations in which a relevant mass change can affect the evolution of the binary system, 
and thus also the waveform of the emitted GW. 

2 Evolution of a binary system in presence of mass accretion 

The evolution of a close binary system in presence of mass loss or gain has been thoroughly 
analyzed in the past, for example as a tool to explain the irregularities and secular changes 
in some spectroscopic binary stars (Kruszewski 1966). Only later a similar analysis has been 
applied to relativistic objects, like a neutron star binary (Clark & Eardley 1977; Jaranowsky 
& Krolak 1992). 

We consider a binary system consisting of two supermassive black holes, following cir- 
cular newtonian orbits around the common baricenter. This system, according to General 
Relativity, radiates gravitational waves, which subtract energy and orbital angular momen- 
tum from the system itself. The effect of this radiation is also to circularize the orbit. For 
this reason we have assumed the orbits to be circular. Also the frictional drag exerted by 
the surrounding gas and stars during the formation and early evolutionary stage of the 
BBHS should circularize the orbits, even if this point is still controversial (Begelman et al. 
1980; Fukushige et al. 1992). 

The angular momentum carried away by the GW in the unit of time is given by the 
quadrupole formula (Peters 1964), which for point masses reads (unit c = G = 1 hereafter) 



where a is the separation between the BHs, Mt = mi + m2 is the total mass, and /x = 
m\m2/MT is the reduced mass of the system. 

When we consider such a system in the core of an AGN, we must also take into account 
the effect of the accretion disk. This material, made up of gas, stars or even more exotic 
objects, like ordinary black holes, is able to transfer or subtract energy and orbital angular 
momentum to/from the BBHS, exerting a frictional drag on the motion of the BHs. We 
will assume that i) the mass falls radially, i.e. without any component tangential to the 
orbits, ii) the direction of the drag force which acts on the BHs is opposite to their velocity 
vectors, and, finally, iii) we neglect any gravitational effect of the disk itself on the BBHS. 

Then a simple newtonian analysis shows that 



For simplicity, we will assume, from now on, that the two BHs have the same mass, m. 
Then, substituting equation (1) in equation (2), one gets the following equation governing 



dJaw _ 32 i?M^ 
dt ~ y a7/2 



(1) 




(2) 
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the separation a 



5 \ a J \m, 



(3) 



In order to integrate equation (3) we should have to know the time dependence of m. 
Unfortunately, while we can make some reasonable hypothesis on the magnitude of m, 
based on the Eddington assumption (see later), very little can be said about the real time 
dependence of rh, as well as the duration of the active phases. Due to our ignorance on 
this subject, we will make the very simple assumption m = = const. Then, integrating 
equation (3) and taking t = as the initial time, one gets 



a{t) = ao IH 



1 + ^ 



where we have introduced the two fundamental time scales 
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16t„ 



(4) 



(5) 



and 



To = 



512 



(6) 



The meaning of these two quantities is as follows: gives the mass e-folding time for 
accretion at the rate u, while Tg gives the time required by the binary to coalesce as a 
consequence of the energy and angular momentum carried away by the GW, assuming 
circular orbits and = 0. For a very slow accretion rate, i.e. for ^ Tg, the time 
evolution of a is dominated primarily by GW emission, and equation (4) reduces in this 
limit to the well known formula (Misner, Thorne, & Wheeler 1973) 



1/4 



a{t) = ao 1 - 



(7) 



From equation (4), we can easily find the coalescing time, i.e. the time at which a 
vanishes. 



'7~r. — Tn 



(8) 



In the limit Tm ^ Tg, by expanding equation (8) in powers of Tg/Tm, to first order we get 

2' 



^ 'Tm. \ 'T'm. 



(9) 
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Figure 1 shows the value of Tc as a function of r^, in units Tg = 1. As expected, the mass 
accretion speeds up the cvohition fo the system. In particular, if the two time scales are 
almost equal, then the coalescing time is approximately 20% of Tg. 

The crucial quantity which appears in equations (4) and (8) is thus the ratio 

^ ^ 3 X lO^mg-^ ( ( , (10) 

Tm ' VI PC J U/©/ yrj ' ^ ^ 

where mg is the initial mass in units of 10*^ Mq. For values of m and m typical of AGN 
models, this ratio is usually bigger than one, unless one is willing to adopt extremely small 
separations. This fact is of considerable importance, because we are forced to abandon the 
usual idealized description, based on equations (6) and (7), in favour of the more general 
equations (4) and (8). On the other hand, the quantity given in equation (10) can not be 
arbitrarily large. In the standard accreting model, the total radiated power is assumed to 
be limited by the Eddington luminosity 

Le ^ 2.6 X 10''^ ms erg s"^ . (11) 

In principle, this limit only applies for spherical accretion, when most of the models actually 
gives L <C (Chang &: Ostriker 1985; Park & Ostriker 1989), but valid arguments suggest 
that this limit remains valid in every realistic (i.e. without an unnatural segregation between 
radiation and fuel, see Turner 1991) anisotropic model, which generally gives L ^ Le (Rees 
1984). 

Associated with is an Eddington accretion rate, that would be able to sustain an 
Eddington luminosity with efficiency e for conversion of mass into radiant energy 

0.2 

ifiE «i — ms Mq/yt. (12) 

Equation (12) gives a lower limit for the e-folding time 

TE ^ 4:A X 10^ e yr, (13) 

which is independent of the mass. Note that equation (12) implies an exponential growth of 
the mass and its derivative, in contradiction with our hypothesis rh = const. Nonetheless, 
since we are implicitly assuming that we observe the AGN for a time T ^ Tm, during 
which the r.h.s. of equation (12) can be assumed constant, we will consider the Eddington 
values (12) and (13) as upper (lower) limits for m and r^, respectively. Thus, the ratio in 
equation (10) is bounded from above by the quantity 



Some comments should refer to the unknown physical parameter e, which describes the 
radiative efficiency. A detailed review of this aspect has been given by Turner (1991), in 
the attempt of accounting for the masses and luminosities of the oldest quasars {z > 4). 
Prom his considerations, one can adopt the reasonable value e = 0.1, independently of the 
details which perturb the fuel reservoir quite far from the inner giant object. 

Finally, another useful quantity is the critical separation, a^, defined as the initial sep- 
aration which gives coalescence after a certain time t. From equation (8), this is given 



In using equation (15), we must remember that the origin of time was taken at the instant 
of formation of the binary. This quantity is of particular importance when one is looking 
for the GW bursts resulting form the final coalescences, since the probability of these events 
depends on the percentage of systems which formed with separation around adtobs)- 

3 The Braking Index 

In the previous section we have calculated the decreasing of the separation a as a conse- 
quence of the emission of GW and the accretion of mass. Since the frequency of the waves 
is twice the orbital angular frequency, the net effect of their emission is an increase of their 
own frequency. If the masses are constant, then a is governed by equation (7), and from 
Kepler's third law one finds 



However, in our case the situation is not so simple, since now the total mass is increasing, and 
a{t) is given by the more complicate equation (4). The overall effect is, thus, a more rapid 
increase of f{t). This fact can have considerable importance in relation to GW detection 
experiments. In fact, a gravitational train emitted by a binary system is detectable as a pure 
sinusoid only if the frequency remains in the same resolution bin during the observations, 
i.e. as long as fT ^ A/ = 1/T, where T is the duration of the experiment. This constraint 
has been sometimes overlooked, in the past, since it implies that the spectral region where 
the signal can be searched for must be consistently restricted (Giampieri &; Tinto 1993). 

In order to describe the time evolution of f{t), we take into account also its derivatives 
/ and /. In particular, the quantity of interest is the braking index 



by 




(15) 




(16) 




(17) 
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Now, Kepler's third law and equation (4) give 



/ 3d Im 3- 21^^/3 5i/ 
7 = -2a + 2™ = — +„• 

Deriving again equation (18) with respect to t, and neglecting the second order terms in 
the nondimensional quantity v{'Kmf)~^/^ , we eventually find 

A; ~ y - /?i/(7rm/)-«/3 , (19) 

where /3 is a numerical coefficient given by /3 2± 0.77. The value k = 11/3, corresponding 
to the limit case 1^ = 0, is the constant value of the braking index when one consider GW 
alone, as one can easily check from equation (16). If wc were able to detect a sinusoidal 
GW at the frequency / emitted by a BBHS, then measuring k one could determine the 
interesting quantity i>m~^/^. Figure 2 shows this result when v is assumed to be equal to 
the Eddington limit (see eq.(12)), with e = 0.1. Given the GW frequency /, or, equivalently, 
the orbital period P = 2//, and the deviation AA; from the general relativistic value 11/3, 
one can determine the mass of the system. Alternatively, if the mass is known, one can get 
relevant information on v, as shown in Figure 3. 

To be more useful, this method ought to be extended to more realistic situations, where, 
for example, mass transfer between the two components occurs. In other words, we would 
like to drop the simplifying assumptions of equal masses and constant accretion rates, and 
find the expression of k which generalizes equation (19). This request can be easily fulfilled, 
since to obtain equation (19) we have not made any use of the solution (4) for a{t); all we 
need is an expression for d in terms of a itself, /x, Mt and their derivatives /i, Mt- From the 
adiabatic invariant one finds 

-2^-#^, (20) 



m 



while equation (1) gives 



64mM| ^^^^ 
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Substituting equations (20) and (21) in equation (18), we obtain, neglecting terms of order 
0{ra/mY and 0{m/m), 



3 96 I /X 3 Mr ' t\jj > \ J 



which reduces to equation (19) when mi = m2 and m = const. From equation (22) we see 
that k, in the general case, can assume values above and below the reference value 11/3, 
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according to the sign of the term between parentheses in the last equation. For example, if 
some material is flowing from one BH (call it component '1') to the other one (component 
'2'), then 



( 



mi — 771,2 



Mr 



(23) 



Mr = 0. 



(24) 



Thus, k is bigger (lower) than 11/3 if and only if mi is less (more) massive than m2- The 
mass transfer has, in practical, no effect on A; as long as mi ~ m^- Moreover, we stress 
the possibility that a single BBHS could exploit different behaviors of ji and Mt during 
its evolution. From the observational point of view, this fact gives rise to a dispersion of 
the measured values of k, depending on which of the various possible effects, namely mass 
loss-transfer-accretion, is dominant on each particular system at that particular time. Thus, 
measuring A; in a wide population of BBHS, one can get very interesting information on the 
variability of these systems, which seems to be an important step toward the understanding 
of the fuel mechanisms. 

4 Evolutionary effects 

Following the discussion concluding the previous section, we will now consider a population 
of BBHS. Independently of their formation epoch and mechanism, we will start considering 
them at a given 'initial' time Iq, when each of them is characterized by a separation oq. 
For simplicity, we will assume that mo and v are equal in all these systems, and focus our 
attention on the evolution of a. In particular, we want to find the distribution function, i.e. 
the function which gives the number of systems with a given separation a, at an observation 
time i, given the same distribution at the initial time to- To be more specific, we define 
the number of BBHS with initial separation between oq and oq -|- da^ as gQ{ao)dao/ a^, 
with ao belonging to the interval [amin, Omax]- The lower limit Omin can be considered as 
the minimum separation compatible with our newtonian assumption (for Post-Newtonian 
corrections to the binary orbit see Lincoln & Will 1990), while the upper limit CLxaax. is 
determined by the sensitivity and the bandwidth of the receiver. 

We want to find how gQ evolves with time. At a given time t, the binary systems with 
ao < ac{t) have disappeared, while those with a > ac{t) are distributed in accordance with 
the number conservation law, i.e. 




(25) 



8 



where a is related to qq by equation (4). Differentiating equation (4) we thus find 



g{a, t) = go{ao) f — ) f 1 + — 



12 



(26) 



GW 



DISK 



In equation (26), we have made exphcit the origin of the multiphcative factors. In fact, the 
term (a/ao)^ is typically due to the emission of GW alone, as can be seen differentiating 
equation (7) (Bond & Carr 1984). Note, however, that in our case the relationship between a 
and ao is not the same as in the 'unperturbed' equation (7). In other words, the factorization 
in equation (26) is only apparent, since appears also, through equation (4), in the term 
indicated with 'GW. Equation (26) can be rewritten as 



g{a,t) = go{ao)F{ao,t) , 
with the introduction of the transfer function F{ao,t), given explicitly by 

t 



F{aQ,t) = 1 + 



16r„ 



1 



1 + 



16 



(27) 



(28) 



Note that this quantity is always less than one, and depends on oq only through the GW 

time scale Tg. Since Tg oc Og, this means that the evolution is much faster near fflmin 
than near Omax- Figure 4 gives an example of this behavior. Note also that, as expected, 
F{acit),t)=0. 



5 Conclusions 

In this work we have considered the effect of an accretion disk on the GW emission from 
a BBHS. The main result is equation (4), which gives the evolution of the separation a as 
due to the GW emission together with a (constant) mass accretion. We have found that the 
coalescing time can be considerably shorter when the mass increases at the Eddington rate. 
Finally, we gave an approximate expression for the braking index, in terms of the reduced 
and total masses of the system. 

We stress that most of the results presented here have been obtained under quite general 
assumptions. Therefore, they can be applied also in other interesting astrophysical objects, 
like a binary pulsar, etc.. For example, we can consider a stellar binary system, loosing its 
mass adiabatically. Then, assuming again equal masses and rh = const. {Tm < in this 
case), we find that a evolves according to the same equation (4). 

In some circumstances, however, the adiabatic hypothesis could appear inadequate. In 
this case, equation (4) is no longer valid. Anyway, we are able to replace it, as long as the 
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fractional change in the binding energy due to the stationary mass loss is proportional to 
the fractional change of the mass itself, i.e. 



SE _ dm 

E m 



(29) 



This case includes the previous one, since the adiabatic law corresponds to a = 5, as well 
as other important cases; for example, Jeans's mode of mass ejection implies a = 3 (Huang 
1963). 

Then equation (3) must be replaced by 



5 \ a 

which can be easily integrated, to give 



128 fm\^ , X , - , 
— (-) -(«-2)(-)m. 



o = ao 1 + 



2-a 



where 



i^a(i) = < 



Correspondingly, the coalescing time (8) becomes 



1 + 



4(a - \)Tg 



t \4("-i)' 

1-(1 + - 

r. 



l_!i!iln(l + 



a 7^ 1 



1 I 4(«-l)r, \ 



. exp(rg/r^) - 1 



l/4(a-l) 



a = l 



When \Tm\ 3> Tg, equation (33) gives 



1 + 



5 — 4a To 



2 r„ 



+ 
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Prom equation (34) we deduce that 

Tc ^ Tg for a ^ 5/4, 

Tc c± Tg for a 2± 5/4, 

Tc ^ Tg for a ^ 5/4. 



(30) 



(31) 



(32) 



(33) 



(34) 
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These formula could be encountered, for example, in a type I Supernova progenitor, 
according to the Double Degenerate White Dwarfs (DDWD) model (Iben & Tutukov 1984). 
In this model, both WD loose their H-rich envelope just before the final coalescence. In the 
usual description of the DDWD model the GW emission is considered only as a mechanism 
capable to get the binary close enough for the common envelope phase (or phases) to 
occur. This idealized model, in which the various phenomena, namely the mass transfer 
and ejection, and the final coalescence, take place with different time scales, is obviously 
justified by the intention of giving a satisfying description of the SNI progenitor. However, 
from the results of our work, we can conclude that all these effects can be taken into account 
simultaneously, and that they can have considerable influence on the observability of this 
kind of sources. 
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Figure Captions 



Figure 1 

The coalescing time, function of the mass e-folding time, r^. Both quantities 

are expressed in units of the GW time scale, Tg. 



Figure 2 

How to determine the mass of the binary, knowing the GW frequency / and the braking 
index k, expressed here in terms of the deviation from the GR value 11/3. We have assumed 
an efficiency e = 0.1. 



Figure 3 

How to determine the efficiency e, for a BH mass of 10^ Mq, from the knowledge of the 
GW frequency / and the braking index k. As in fig.2, AA; = 11/3 — /c. 



Figure 4 

The evolution of the distribution function giao) at a generic time t. Binaries with 

do < o-c{t) have disappeared. Those with an initial separation larger than adt) have evolved 
according to equation (27). The initial distribution go is assumed to be flat over the interval 
[0, Omax] (dot-dashed line) . 
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